Table of areas for bounding box sites (area in square Hectometers)

Histogram of areas in square meters, separated in to size classes to make easier to visualize.

Below is a table of all the bounding box sites on (or intersecting) federal Indigenous lands.

## Spherical geometry (s2) switched off
## although coordinates are longitude/latitude, st_intersects assumes that they
## are planar

Are sites that intersect reservations differently sized from sites that don’t?

## 
##  Welch Two Sample t-test
## 
## data:  as.vector(on_rez$area) and as.vector(not_on_rez$area)
## t = 6.213, df = 231.08, p-value = 2.398e-09
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  13287023477 25627840106
## sample estimates:
##   mean of x   mean of y 
## 19587555177   130123386

The graph below summarizes areas by datasettypes, only for those datasets for which there are more than ten bounding box sites.

## although coordinates are longitude/latitude, st_union assumes that they are
## planar

Dataset types on federal Indigenous lands versus others?

## Joining with `by = join_by(site_name, siteid, area, size)`

## although coordinates are longitude/latitude, st_union assumes that they are
## planar

## corrplot 0.92 loaded

## `summarise()` has grouped output by 'datasettype'. You can override using the
## `.groups` argument.
## although coordinates are longitude/latitude, st_union assumes that they are
## planar